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WRINKLING PHEFOMEICA 'OF THIN FLAT PLATES 
SUBJECTED TC SHEAR STRESSES* 
By F. Bollenrath 

The present report covers a series of tests on thin, flat 
elastic strips restrained at two parallel edges and subjected 
to shear by conversely directed stresses. The stress is uni- 
formly distributed over the constrained edges and acts in the 
plane of the strip. As soon as the shear in such strips exceeds 
a certain limit the equilibrium between the applied and .the in- 
ternal stresses is interrupted and instability prevails. Now, 
if the wall is thin with respect to the width of the median 
plane of the body, the incipient instability is attended by a 
characteristic flexural phenomenon which aopears as a pronounced 
lobe or wrinkle. Thus, Figure 10 shows two restrained duralumin 
plates where the type of shear just mentioned, resulted in a 
wavelike lobe. The deflection in these plates is so pronounced 
as to remain visible when the load had been relieved. 

Knowledge of the conditions under which this wrinkling under 

shear (i.e., the shear producing instability) occurs, is of great 

import for many structural components, particularly in light 

metal structures. The problem treated here is applicable to all 

cases where a thin wall, through the type of its attachment 

* 11 Ausb'euler scheinungen an ebenen, auf Schub beaaspruChten Flat- 
ten.* From Luftfahrtforsolxung, Vol. 71, Ho. 1, December 12, 
1929. 
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to other pa.rts, is intended to take up and transmit stresses, 
as in plate girders, bulkheads, fuselages and wings of various 
airplane types. 

Two limiting cases (Fig* 1) for thin bodies stressed in 
shear are those of round plates and thin-walled tubing (Refer- 
ence 17), the latter of which, however, really belongs to the 
"shell * group. The first experiments on this subject were 
made by Lord Rayleigh (Reference 14), who investigated the 
case of a vibrating system with reference to its natural fre- 
quency of vibration. Bryan (Reference 8), Walter Ritz (Refer- 
ence 15), and S. Timoschenko (References 21 and 22) based their 
experiments on the stability equations of an ela.stic plate for 
defining its buckling stresses.* 

In the following we describe the development of the calcu- 
lation of this problem from approximation to exact solution. 

Let E = Young T s modulus (kg/cm 2 )', 

m = transverse elongation factor (Poissoi^s ratio), 
M. - 1/m, 

L = length of rectangular plate (cm) 

h = thickness of plate (cm), 

• a = width of plate (cm) , 

U| v, w = displacements in the direction of 
x, y, and z axes 

*St. Bergmann and H. Reissner reported, on the appearance of in- 
stability in corrugated plates under shear in n I T euere Probleme 
aus der Flugzougstatik, " end ,f uber die Knickung von Wellblech- 
streifen bei Schubbeanspruchung. ,r Zeitschrift fur Flugtechnik 
und Liotorluftschiffahrt, No. 18, 1929, pp. 475-481. 
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shear per unit of length in restrained 
edge (kg/om), 

compression per unit of edge surface parallel 
to x axis (kg/cm 5 ) , 

tension per unit of edge surface parallel to 
y ar.is (kg/om 2 ), 

p/q> 

Euler's Duckling load for a strip of length a, 
thickness h, and width 1 (kg/cm 2 ), 

wave length parallel to the restrained edges, 

1 ... flexural rigidity, 

12 (1-M- 3 ) 

a 5 s 
D 

Wrinkling Due to Shear (W. S. Lilly's method) 

Lilly (References 10, 11, and 12) in his work on wrinkling 
of plate girders, applied the well-known strip method. Thus, 
we see in Figure 2, the wavelike deflection of the median plane 
of a plate under shear after exceeding the shear strength, ex- 
pressed as lines of equal deflection perpendicular to the me- 
dian plane. The figures attached to the contours have a rela- 
tive significance only. We car. visualize the same deforma/tion 
as having been induced "by tension and compression stresses of 
the kind shown in Figure 2, and which are dimensioned in magni- 
tude and direction so as to balance a shear stress of the in- 
tended type in magnitude and direction. The applicability of 
Lilly's method of approximation rests on this theorem. He con- 



s = 

P = 

q = 

r = 

o e = 

X = 

D = 

A = 
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eiders, according to Figure 3, one single panel - between two 
stiffeners - of a web member collapsed in the form of a half- 
wave. The line o-o denotes- the center line of the half-wave. 
Such panel then is stressed in tension and compression as a con- 
strained rectangular plate of length L and. width X/3 (See 
Fig. 4)./ r 

How when we divide the plate in strips b-b perpendicular 
to 'each other, parallel to axis at. and a~a, and parallel to 
axis y, strip a-a is stressed in bending and compression, 
and strip b-b in tension and bending. So when we assume a 
sinuous deformation and apply Euler ? s formulas for the buckling 
load we find 

P ^ r L 2 / 12 LV' UJ 

There is a minimum for equation (l) when 



^ - J ^ + 1 = J . (2) 



4 If r 

This formula (2) shows that the wave length is not affected by 
the plate width; it merely depends on L end r. Formula (l) in 
connection with the \ values of formula (2) yields 

571? " / islf + i -~T 

If, by L = \/2, the wave runs at 45° to the constrained 
edges/ the stress is a pure shear stress for r = 1. In that 
case, formula (2) yields 
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\ 2 = A {fz- 1) L 2 - 1.64 L 2 



(4) 



According to these formulas the applicability of such calcu- 
lations is contingent on the condition that the dimensions of 
the plate allow the formation of a half-wave only. 



This approximation method is based on the potential energy 
of the system (Reference 21). The appearance of the unstable 
equilibrium attitude becomes evident as a potential energy in 
bending V, when deflecting from the plane, and as a lower po- 
tential energy In shear \ , so that 



Starting with a rectangular, all-around supported plate, 
the length-width ratio L/a is progressively increased and soon 
reaches a limit value, valid for the infinitely long plate strip 
the type of support of the smaller sides becomes secondary. In 
this manner the formula for the wrinkling of a long plate from 
the median plane, w becomes: 



Wrinkling Due to Shear (S. Timoschenko ' s method) 

(Reference 22) 



(5) 



= A sin tt £ sin tt 

D 



| (x - a y) 



(6) 



and for wrinkling due to shear 



t = £e [ 6 a 2 + 2 + — ^- + kj£ (l + a= f ] 




(7) 
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For a minimum' t we Have : 



X 
2 



a 



J'l + a 2 



(8) 



a 



1 



(9) 




so that 



T 



(10) 



Figure 5 shows the c values for various relations of 
plate length L to plate width a. It will he noted that c 
soon approaches the value for L/a = oo . 

Wrinkling. Due to Shear (R . V. Southwell end Sylvia W. Sfcan) 

These authors give a strict solution (References 18 and 19) 
for the case of a flat elastic strip of infinite length. But 
the distribution of stresses in a thin plate over an infinitely 
small plate element is without any pr:.?tiaal significance, ac- 
cording to Love. For that reason we represent the stresses ap- 
plying at the edges by their resultants and their resultant mo- 
ments with respect to the unit length of the border line, i.e., 
the curve in which the edge passes through the median plane. 
How according to St. Tenant's theorem the effects of two systems 
of; stresses yielding.. .identical resultants and resultant moments 
are practically identical at some distance from the edge. 

Selecting from Figure 6 a plate element dx long and dy 
wide, we have,. .since the stress in a plate depends on the shape 
of the elastic surfa.ce: 
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1. On the edges parallel to the x axis - To ending stresses 
linearly increasing with the distance from the median plane, 
statically equivalent to a force couple with moment a, ; 
t orsional stresses equivalent to a torsion moment Hj^ J shear 
st resses in the direction of the z axis with resultant H s J 
shear st resses in the direction of the x axis with resultant S. 

3. On the edges parallel to the y axis - bending stresses 
linearly increasing with the distance from the -median plane, 
statically equivalent to a couple of moment Cr 1 ; torsion 
stresses equivalent to a torsion moment H^J sh ear stresses in 
the direction of the z axis with resultant N x J s hear st resses 
in the direction of the y axis with resultant S. 

How the equilibrium equations applied to the deformed plate 
in Figure 7 are as follows: 



a * ay 




- o 



(ii) 



ax ay 
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d 2 




9 l\ 


9 x 2 


By 3 x 










+ 


a £ 


9x By 




ay 


0 1*1 






a 



= 0 



(12) 



9 3x ay 



2 a 

^2 

9 Y 2 



9" W 
+ 3 s 

ax ay 



c 



j 



The stress resultants and the result ant stress moments are 
combined to read: 



v a x 2 ay 2 / 



^9 y a 9 x 2 / f 



H x '« D (1 - M-) i 



P* 2 TT 



9x ay 



(13) 



This, together with the above formulas, yields 



:: (¥JL + u La£jL^ + 2 (i - W -Jf 



w 



v -3jc 4 



9x 2 9y 2 



+ 



o* 9y' 



a y 



3 _ 2 / 



ax ay 



V 9x 4 



9 4 w + a* w\ 



9x 2 9y 2 9 y 4/ 



9 2 w 



3 S 



3x 9y 

a 8 w 
ax a y 



D A A w + 2 S ^ " = 0 

^x ay 



(14) 

(15) 
(15) 



Southwell resolved this last equation (16). Writing 



w = i e 



, (i - 1) 



(17) 
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formula (16) becomes 

v(*) ~ 3 k 2 y ,! + 2 i * S Y ! + * 4 Y = 0 (18) 

for function Y (y) , and. must be integrated for the limiting 
equations 

y -e ± | r Y - Y« =■ 0. (19) 

Finally n v . 

Y = e lA,y (20) 

yields for X an equation of the fourth degree, Whose four 
roots give the solution. According to this solution the median 
plane of the plate deflects, due to exceeding the shear strength, 
into a uniformly loted surface. The lowest shear strength is 

s m 88.7 % kg/cm (31) 
a 

and the conformal wave length. 

\ * 1.6 a (32) 

This length is contingent on the plate width, as shown by 
formula (2). Thus Figure 8 presents the wrinkles due to shear 
plotted against the ratio: wave length/plate width, according 
to Southwell and Skan." The curve of the first solution is ex- 
pedient for the case of deforming into one, the curve of the 
second solution into two rows of half-waves in the direction of 
the plate width. Accordingly, a plate under simple shear, whose 
stiffeners bring about two rows of half-waves in the direction 
of the plate width, is more conducive to higher shear strength 
than one producing but one row. 
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S x p e r i m e n t s 

To produce a sheep stress in a plate strip of the kind de- 
scribed, we developed a clamping arrangement which constrains 
the test strip on two parallel edges at a width of 5 mm between 
flat surfaces, as exhibited in Figure 9, where (a, a,) are the 
test strips and (b,b) the two constraining clamps. The 
stresses acting along the constrained edges shifted the latter 
parallel to each other and put the strips under shear. 

The loading was accomplished by applying a tension on the 
inside clamps, which is evenly distributed over both test strips 
and taken up by the outside clamps. 

For the experiments on brass and duralumin we used a 5- ton 
tension machine (Fig. 10, vertical arrangement), and for the 
celluloid tests we adhered to the horizontal method (Coker) be- 
cause of the small loads needed and the comparatively heavy 
weight of the clamping device* 

The relative displacement of the constrained edges was 
measured with four Zeiss gauges (one division = l/lOO mm). 

Conditions of the Boundaries or Edges 

The problem is whether and to what extent the conditions 
and assumptions upon which the calculations are based are com- 
plied with in these tests. Southwell and 6k an analyze, for ex- 
ample, an infinitely long plate strip, a condition which is prac- 
tically unrealizable. The strips used here are of limited length 
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and with free edges; consequently the problem is to discover 
what effect this substitution of finite- strip length for one of 
infinite length' relative to the stress distribution over the 
whole strip has on the magnitude of the shear stress and on the 
type of deformation. Then, another question to be answered is, 
the role played by the type of support of the smaller sides or 
of the ends of the strip. 

Professor E.G. Coker (Reference 5) reports on his optical 
experiments with transparent bodies subjected to shear. He de- 
fines the stress distribution along a line '0 X (Fig. 13). The 
boundary conditions for the plates are: The long sides are 
constrained and are shifted parallel in the plane of the plate; 
the short sides are free. Some of his results ere indicated 
in Figure 13. The origin of the coordinates is in the center 
of 0 X. The plate lengths (in inches) are the abscissas and 
the stresses (in pounds per square-inch) are the ordinates. 
The data in Figure 13 are for a plate 2 inches wide. According 
to his tests the shear distribution in a rectangular plate is 
parabolic along 0 X (curve a in Fig. 13). When the ratio of 
a to L decreased the curve of the stress distribution first 
flattened out at the tip (curves b and c, Fig. 13). For still 
smaller A/L ratios, local raises in stress distribution oc- 
curred near the free edges (point 1 on curve d and point 2 on 
curve e). • In the middle of the strip the distribution w a s 
practically uniform. 
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Anclrade (Reference l) gives an experimental and analytical 
treatise of the problem. His tests were made on a lb ,r x 4" x 4 M 
rectangle of gel, a glycerine-gelatine compound. The two long 
sides were nailed by thin copper screening to a "board. Then the 
screening was heated and pressed on the gel-body. The deforma- 
tion was produced by parallel shifting of the plates. For de- 
fining the displacement he applied the optical method of Charles 
Pearson and A. F. 0. Pollard ( !l An Experimental Study of the 
Stresses in Masonry Dams , Drapers Company Research Memoirs , Tech- 
nical. Series II, V) . Andrade observed the same local stress in- 
creases near the free edge, but his calculations do not agree 
with his experiments, a fact which prompted him to speak of the 
i na.de qua cy of his applied analytical methods. 

L. N. G. Filon (Reference 6) examined a beam of rectangular 
cross section which is stressed in shear applied to the upper 
and lower fibers. Here 0ol:er T s observations are strikingly 
confirmed. 

S. Inglis (Reference 7) treats the same problem; but attacks 
it from a different angle and his results axe in close agree- 
ment with Professor Ooker's i'ijures. Inglis based his calcula- 
tion on the superposition of two cases. In the first case he 
constrains the plate (Fig. 14) at ed^es A3 and CD, leaving 
AD and BO free. Then, he subjects the edges to a uniform shear 
acting in. the plane of the plate. In the other case, the plate 
AB,0D is so stressed in shear, which is uniformly distributed 
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over the four edges , 'that the 'shear at the edges AD and BO is 
in opposite direction to the shear of the first case. The re- 
sult is an absence of shear on edges AD and BO. 

The results of •Inglis' experiments are exhibited in Figures 
14-18. Figure 14, in particular, shows the distribution of 
shear along OX (y = o); the edges AB and 0D are constrained, 
the shear is uniform over AD. In Figure 15, we note the shear 
perpendicular to boundary A3 (y = 1); it is uniformly distrib- 
uted over AD, with AB -and 0D constrained. Figure 16 shows 
the almost parabolic distribution along OX (y * o) for a 
square plate; Figure 17, the shear along OX for a plate with 
ratio length - width = 2 ; and Figure 18, for ratio - length/ 
width = 4. In the last figure the flatness of the shear curve 
is very noticeable. The effect of the free edge doec not last 
beyond 1.5 times the plate width to plate center. In conformity 
with the tests of the author, the smallest ratio of plate length 
to plate width = 8.5, justifies our assumption that the stress 
distribution is practically uniform over the whole strip, and 
that the free edges have no material effect on the test data. 

The experiments included tests with celluloid, duralumin 
and brass. 

Celluloid 

For the calculation of the deformations and of the wrinkles 
due to shear, it is absolutely essential to know the longitudi- 
nal and the transverse elongation factors. For this reason, 
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special stress was laid on the determination of these figures, 
which was accomplished by micrometric measurement of rectangular 
test plates with the Martens reflecting instruments. The length 
of the spring was 100 mm for the longitudinal, and 25 mm for the 
transverse elongation. By a distance of 1135 mm and a height 
of 4.5 ram, a 1/5000 gear ratio was obtained. The tests were 
made on the 10-ton Losenhausen machine in the Aerodynamic Insti- 
tute of the technical high school at Aachen, which also included 
the use of the 1000-kilogram metering box. 

Celluloid proved particularly- suited to such experiments. 
The proportional elastic limit was reached at 130 kg/cm 2 tension 
and by a comparatively high total elongation of 0.6^. At fail- 
ure the material showed a 7-10 per cent elasticity by a total 
elongation of 26 per cent. 

We examined eight different kinds of celluloid, according 
to Table I. The stres-3 strain curves in Figure 19 show the 
characteristic behavior of the celluloid until failure: an ap- 
proximately even raise in elongation with the stresses up to 
the yield limit, above which the breakdown becomes very percep- 
tible. There is a basic difference between the clear, trans- 
parent and the colored celluloid beyond the yield point. The 
colored material, shows a marked drop in stress, followed by a 
subsequent raise. This is most likely due to the coloring mat- 
ter . 

The nature of celluloid to follow Eooke T s law to the yield 
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point and then to enter a so-called plastic state, proved very 
fortunate in subsequent wrinkling tests. In these we used the 
plates 5-8 of Table I, which, with respect to modulus of elas- 
ticity (S == 15,800 to 23,300 kg/cm 2 ) and the transverse elonga- 
tion (m =? 2.415 to 2.8) gave us more room for action. 



TABLE I. Modulus of Elasticity of different Kinds of Celluloid 







E 










Modulus 


Transverse 




Number 


Kind 


of 


elongation 


Thickness 






elasticity 


h 






kg-/cin 2 


m 


cm 


1 


clear 


36,600 


3.4 


0.19 


2 


fi 


' 35,700 


2.4 


0.19 


3 


white 


24,500 


3.0 


1. 51 


4 


n 


23,500 


3.66 


0.08 


5 


clear 


28,000 


3.45 


0.053 


6 


fr 


16.000 


3.8 


0.05 


7 


if 


23 , 500 


3.56 


0.O5 


8 


white • 


33,300 


3.45 


0.031 


9 


clear 


13,500 


3.415 


0.051 



Particular attention was paid to temperature and humidity while 



making the tests. 

Test Data 

a ) Wr i nk 1 i n g s t r ess.- 3y a certain shear the constrained 
edges of a plate width a may be shifted parallel by an amount 
g 1 . Expressing a and b in centimeters, formula 

g = Si 100$ (23) 
a 

now is to denote the shifting of the plate width in per cent. 

The behavior of the plate subjected to shear is analyzed 
in Figure 20, where the loading s (kg/cm is plotted against 
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the raaacimum deflection from plan© w (in centimeters .10 ) 

and against displa.cer.ient g (in per cent). The dimensions of 

• a = . 

the plate were: thickness, h = 0.05 cm; <;.idth,/l.9 cm; and 
length, L - 23.7. By critical shifting here when g = 0.8$ 
(point A), the load deflection line begins to deflect, the load 
increases progress more slowly than before by corresponding dis- 
placements. This incipient deflection of the load line is si- 
multaneously -accompanied by a deflection from the plane. But 
even before the deflection from the central plane assumed meas- 
urable significance, the incipient buckling became noticeable 
through. the distorted image of the reflected plate. 

A second test with a celluloid strip of thickness h =0.05 
cm, and width a = 1.6 cm, is illustrated on Figure 21. The 
load deflection line is straight to point A, which corresponds 
to a displacement g = 0.78$, where a bend occurs. This means 
the proportional limit has been reached and exceeded. Above 
this bend the line continues almost linearly, except for the 
more pronounced slope toward the axis of displacement. But 
the second portion itself shows at B a sudden, stronger deflec- 
tion toward the axis of displacement; the load stages for ob- 
taining equivalent displacements as before, become materially 
smaller. Bulging sets in; the wrinkling shear has been reached. 
But it occurs in this case above the yield limit; the deflec- 
tion is inelastic in contrast to Figure 20. 
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TABLE II. Tests 1-9 

Material , Celluloid 

Modulus of elasticity, E = 16,000 kg/cnf 2 

Poisson 1 s ratio, m = 2.8 

Plate thickness, h = 0.05 cm 



Test 
No . 


Plate 
: width 
a 

ora 


a/h 


wrink- 
ling 
load 

kg/ cm 


Displace- 
ment 




k' 

mean ■ 


k 


Wave 
length 
X 

cm 


k 
a 


1 


1.9 


38 


2.635 


0.600 


76,900 


76,600 


50.0 


4.16 


2 


.18 


2 


1.9 


38 


2.635 


0 . 643 


76,900 




4.20 


2 


• 18 


3 


1.8 


36 


2.945 


0.667 


76,300 


Elas- 




3.68 


2 


.04 


4 


1.8 


36 


2.945 


0.667 


76,300 


tic 




3 .69 


2 


.05 


5 


1.7 


34 


3.31 


0.824 


73,500 


def or- 




3.25 


1 


.9 


6 


1.7 


34 


3.31 


6 . 848 


76,500 


roation 




3.46 


9 

to 


.03 


7 


1.6 


32 


3.49 


0.900 


73,200 






3 .104 


1 


.94 


8 


1.5 


30 


3.62 


0.734 


66,800 






2.91 


1 


.94 


9 


1.4 


28 


3.86 


0.843 


62,000 






2.63 


1 


.88 



By observing the appearance of deflection above as well as 



below the limit of elastic displacements, we made a series of 
tests with celluloid and defined the wrinkling load (Test series 
Nos. 1-9, Table II) . 

We considered as wrinkling load that load at which the 
first sign of measurable deflection from the central plane was 
noted. Hereby it should be taken into consideration that the 
load at wrinkling will be about as high as point B in Figure 20, 
but even at that it- will be 25fj lower than the theoretical figure. 

Thus the quality of the approximation to the theoretical 
buckling load remains to a certain extent a- question of the con- 
straint. The experimental values should prove useful for the 
designer, because they supply a practical criterion of incipient 
wrinkling, for the danger . of collapse lies not so much in the 
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formation of wrinkles as in the material displacements and in 
the pronounced uneven plate stresses accompanying the inception 
of deflection. 

The results of the first test series (tests 1-9) are given 
in Table II and Figure 22. The curve representing the depend- 
ence of the wrinkling load on 1 a/h may be expressed in the 
form of the equation 

s = k 4 (33) 
a 

When a/h falls below a certain value the experimental 
wrinkling loads with decreasing a/h deviate more and more from 
those expected according -to the preceding formula. The basis of 
this limited applicability is Euler T s load formula* Back in 
1845 E. Lamarle (Reference 9) pointed out that the proportional 
limit of Euler ! s formula in its original form was to be accepted 
as limit of validity. (Todhunter and Pearson, "A History of 
Elasticity and Strength of Materials , u Vol. I, No. 1253 ff 
Cambridge, 1886.) ' The unlimited validity of Euler ! s formula in 
a more general form was established by Engosser (Reference 3), 
and appeared' in 1395 (Reference 4), the accuracy of which was 
later confirmed by v. Ear man (Reference 8). 

The proportional limit is exceeded whenever wrinkling occurs 

under a smaller load than could be expected according to formula 

s = k -2-. The correct relation between shifting and stress, upon 
a 2 

which the calculation is based, exists no longer. As the a/h 
ratio becomes smaller the assumption of a "thin plate" loses in 
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identity. Expressing formula (23) 



as 




1 




|JT 12 



k ! remains, as long as the deflection occurs below the propor- 
tional limit for shear, an Invariable as proved "by tests 1-6. 
If, to arrivd at WS inkling, the proportional limit is exceeded, 
k ! becomes so much smaller as the permanent deformation is 
stronger at buckling. : Consequently, k ! yields k from the 
experimental values vrithin the elastic range', vrhich, according 
to tests 1-6, is k = 50. 

Three other series of tests on celluloid, Hos. 8, 9, and 
10, Table I) - tests Bos. 10 to IS, 17 to 26, and 27 to 36 
yield practically the same data, and are exhibited in Tables. 
Ill, IV, and V, and Figure 22. k - 51. 52, and 49.3. The tran- 
sition from elastic to inelastic deflection is characterised by 
a/h - 34, 32.2 to 35.5 and 31.4 to 29.4. 
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TABLE III. 



Material , 
Modulus of elasticity, 
Poisson T s ratio, 
Plate 



Tests 10-16 
Celluloid 



E 
m 
h 



22,500 
2. 56 
0.05 0] 



kg/em 5 



Tost 
llo . 


Plate 
width 
a 

cm 


a/h 


Wrink- 
ling 
load 
s 

kg/cm 2 


Displace- 

ner t 

' $ 


k' 


k 1 

, mean 


k 


10 


1.9 


38 


3.94 


0.82 


115,000 


111,402 


51 


11 


1.9 


38 


3.73 


0.73 


109,500 






12 


1.8 


36 


4.17 • 


0.956 


110,300 






13 


1.7 


34 


4.82 


1.0S 


110,300 






14 


1.5 


30 


5.17 


1.235 


95,400 


Inelas- 




15 


1.6 


32 


4.83 


C.59 


101,400 


tic 




16 


1.6 


32 


5.04 


1.09 


105,800 
















1 


1 illation 





Wave 
length 

\ 

cm 



3.62 
3.62 
3. 38 
2.75 
3.10 
3.01 



TA3LS IV. Tests 17-26 



Material , 

Modulus of elasticity, 
Foisson 1 s ratio, 
;e thickness, 



Celluloid 

S = 23,30C kg/cn J 

m =2.45 

a = 0-031 cm 





Plate 




Wrink- 








width. 


a/n 


ling 


Displace- 


k ' 


Test 


a 


load 


ment 


llo. 






s 








cm 




kg/cm 2 






17 


l.S 


61. 3 


i.bi 


0.43 


121,300" 


18 


1.8 


58.1 


1.075 


0.406 


116,900 


19 


1.7 


54.9 


1.235 


0 . 365 


119,900 


30 


1.5 


48.4 


1.57 


0.502 


118,900 


21 


1.4 


45.2 


1.36 


0.642 


113,400 


pp 


1.3 


42 . 0 


2.14 


0.3 


121,600 


23 


1.2 


38.7 


.3.42 


0.918 


117,200 


24 


1.1 


35.5 


2.95 


1.05 


119,700 


25 


1.0 


32.2 


3 . 36 


1.36 


112,300 


36 


0.8 


25.8 
_ .. 


3.73 


1.81 


72,100 





mean 
1197200 



Inelas- 
tic 
defor- 
mation 



V 


Wave 
length. 

X 

cm 


X 
a 


to 


3.73 


1.96 




3.53 


1.96 




3.16 


1.86 




2.97 


1.97 




2.59 


1.85 




2.45 


1.80 




2.33 


1.93 




3.22 


2.00 




2.15 


2.15 




1.61 


2.01 
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TABLE V. Tests 27-56 

Material, Celluloid 

Modulus of elasticity, E = 15,300 kg/cm 2 

Poisson ! s ratio, m = 3.4-15 

Plate thickness, h = 0.051 cm 



Test 
No. 


Plato 
width 
a 

cm 


a/h 


Wrink- 
ling 
load 
s 

kg/cm 


Dis- 
place- 
ment 

i ' 




k« 

mean . 


k 


Wave 
length 
\ 

cm 


\ 
a 


27 


1.25 


24.6 


5.85 


1.595 


68,900 


in 






2.44 


1.955 


28 


1.3 


25.5 


5.34 


1 . 295 


68,000 


elas- 






2.52 


1.94 


29 


1.35 


23.6 


5.08 


1.125 


69,300 


tic 






2.66 


1.965 


30 


1.4 


27.45 


4.83 


1.13 


71,400 


defor- 






2.76 


1.97 


31 


1.45 


28.45 


4.32 


0.98 


68,500 


ma- 






2.78 


1.92 


32 


1.5 


29.4 


4.32 


1.00 


73,300 


tion 






2.93 


1.95 


33 


1.6 


31.4 


4.G7 


0.99 


79,300 








3.17 


1.98 


34 


1.7 


33 .3 


3.61 


0.847 


80,400 


78,300 


4-9. 


n 


3.38 


1.99 


35 


1.8 


3-5.3 


3.16 


0.722 


77,300 








5.43 


1.94 


35 


1.9 


37.25 


2.80 


0.60 


76,200 








3.8 


2.00 



To) Type of deformation.- Aside from the load which induces 

wrinkling, the type of deflection is of particular interest 



with respect to the wave length, for it enables us to decide 
where reinforcements are necessary. 

We used a Zeiss instrument with which we practically cali- 
brated the deformed plate. 

Owing to the aforementioned irregular stress distribution 
near the free edges, we naturally expected some disturbance in 
the uniformity of the deflection. For that reason we had to 
ascertain the extent of the effect of the free edges on the 
deformation, because this boundary zone had to be eliminated 
before defining the wave length. In several specimens the 
first signs of wrinkles appeared at 1 to 1.5 a from the free 
edge. But by a very slight load increase the deformation spread 
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evenly over the entire strip. . This undoubtedly is due to the 
slightly higher stress increase near the edges. 

So in order to define this edge effect, we examined the 
■boundary zone of a transparent strip of celluloid 0.05 cm thick 
and 1.9 cm wide after wr ink,ling. Our method consisted in deter- 
mining the deformation as dissections perpendicular to the' con- 
strained edges. The dissections were 2 mm apart. The measuring 
points w.ere 1 mm apart. > . : " . '. ' • 

Since the shape of the lever made the measurement- of the 
deformation on the constrained edge itself- impossible, the ends 
of the dissections had to be determined afterward. The deforma- 
tion over several wave lengths was measured 1.5 mm from the con- 
strained edge on a line parallel to this edge. The middle lines 
between the tangents to the waves then yielded the position of 
the constrained edges, after which the contour lines, i.e., the 
lines of equal deflection from the plane, give us the deforma- 
tion, ' shown on Figure 23. 

The figures on the contour denote the deflection perpendic- 
ular to the median plane in 0.01 mm. The fourth half-wave al- 
ready has exactly the same shape as the one following. The 
greatest deflection on the whole strip occurs on the free edge. 
The irregularities in the deformation do not extend beyond 
three half-waves. 

So it may be stated that the effect of the free edge on the 
deformation does not exceed 1.5 times the width of the plate 
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from the free edge, and all waves outside of this defined "bounda- 
ry zone are included when defining the wave length. 

However, since this does not preclude the possibility of 
any more extended edge effect under increasing deflection we 
determined the location of the zero deflection lines or of the 
junction lines for s « 4.83, 5.43, and 5.93 kg/cm at vari- 
ous degrees of deflection. The strip was 1.8 cm wide and 0.05 
cm thick. But the test failed to show 'any difference in dis- 
tance between two junction lines for these load stages, with 
the exception of their slope which, owing to the greater dis- 
placement of the constrained edges, was more pronounced (Fig. 24). 

To follow the course of the wrinkles, we determined the 
deformation on the sajne plate for a half-wave, the results of 
which are plotted in Figures 35-27. At the left we find the 
deflection from the plane as dissection perpendicular to the 
constrained edges and the absolute values of the deflection; 
at the right, the deformation as contour. For s = 4.83 kg/cm 
the deflection is perfectly elastic, while at higher loads part 
of the deflection is already of a permanent nature. From these 
deformations the disproportionate raise in wrinkling by small 
load increases, becomes apparent. The permanent deflection 
after releasing the load is shown cis dissection and contour on 
Figure 28. 

The wave length was defined by measuring the distance of 
maximum deflections from each other. The results of these 
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measurements a:re given In Tables II to V, w^tK'wave length X in 
centimeters . 

, For X/a we obtained the following figures: 

Tests 1 - 9 r ~" I- 88 - 2.18 

10 - 15 1.31; - 3.00 

17 - 36 1.35 - 3.15 

27 ~ 36 1.92 - .2.00 

There is no difference in wave length for elastic or inelastic 
deflection as far as we could ascertain. The average X/a = 
1.S7, as seen in Figure 29, shows a plotted against X/a. 

o ) Stress dist ributio n at t he beginning of wrinkling . ~ Her e 
we attempted to determine the bending moments perpendicular and 
parallel to the- constrained edges from the deflections developed 
in the plate after exceeding the wrinkling load. The plate was 
0.05 cm thick, 1.3 cm wide, of celluloid No. 7, Table I; the 
measurements were made by s - 4.83 kg/crn edge loading. 

Our procedure was to approximate the dissections parallel 
arid perpendicular to the constrained edges as arbitrary functions 
by a series of other given functions. Inasmuch as these dissec- 
tions are valid as periodic functions, we chose "file approximation 
by Fourier series. Then we transferred the measured values to a 
system of coordinates and plotted a smooth curve which followed 
the experimental points as closely as possible. This was fol- 
lowed by an harmonic analysis of empirical functions, after 
which we took a number of equidistant values from the curve. 
The individual periods were divided into twelve parts and the 
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coefficients of the sine and cosine terms of the Fourier series 
computed. Terms with coefficients, which, in their order of mag- 
nitude west below the measuring accuracy of the Zeiss instru- 
ments, were disregarded in the second differentiation. From 
the quotient of the latter we determined the respective 12 ordi- 

nates which yield the values for and that is, for 

ox oy c 

the moments. 

Table VI contains the perfectly elastic deflections perpen- 
dicular to the constrained edges. 

T^e distribution of the computed moments 



over a half-wave may be seen in Figures 30 and 31. The lines ap 
ply to Q % and G 2 = constant. The ascribed figures yield with 
- 1.64 x 10~ 3 as scale the moments G x and G 2 (jcg/om) and 
with - 3.92 the stress in the plate in kg/cm 2 . The maximum 
moments and stresses are in the center of the half-wave and 
amount to 0.469 and 0.335 kg/cm, and 112 and 78 kg/cm 2 . 

It is seen that the plate stresses, due to wrinkling in 
bending, exceed the shear stresses of 96.6 kg/cm 2 at the edge 
soon after exceeding the buckling load. Because of the very 
pronounced increase in bending stress by further deflection, it 
becomes evident that a stress is quickly developed by the wrink 
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les above the proportional elastic limit. 

TABLE VI 

Deflections in Dissections Perpendicular to the 



Cons trained Edges in l/lOOO mm 





a 5 




d 


e7 




* 5 


3 


"k 5 


m 




P 


1 


-12 


-34 


-66 


-93 


-117 


-125 


-128 


-113 


- 87 


-49 


-33 


2 


n 

— U 


-12 


-29 


-51 


- 76 


-105 


-125 


-130 


-116 


-79 


-35 


3 


20 


23 


27 


21 


4 


- 37 


- 68 


- 95 


-102 ' 


-76 


-36 


4 


47 


77 


90 


90 


33 


. 53 


15 


- 37 


- 61. 


-60 


-36 


5 


57 


105 


13o 


149 


148 


131 


82 


3 7 


- 3 


-27 


-34 


6 


60 


110 


160 


179 


301 


197 


159 


105 


46 


13 


- 4 


7 


50 


84 


157 


205 


337 


0">Q 
cj i-' 


216 


167 


100 


43 


15 


8 


35 


73 


129 


182 


2 2 5 


238 


235 


197 


137 


74 


39 


9 


o 
o 


33 


75 


127 


169 


193 


204 


193 


140 


100 


39 


10 


-11 


-10 


+35 


74 


119 


163 


186 


134 


155 


105 


45 


il- 


-38 


-45 


-38 


- 5 


36 


75 


104 


113 


108 


76 


35 


ia 


-49 


-65 


-72 


-67 


-38 


0 


30 


52 


57 


45 


33 


13 


-53 


-73 


-87 


-96 


-93 


- 70 


- 30 


7 

i 


21 


35 


13 



Tests with Duralumin 



These experiments should be of special interest to the user 
of light metal. We used rolled plates 0.05 and 0.03 cm thick, 
alloy No. 661b, hardness 1/2. The elastic properties, modulus 
of elasticity, and Poisson'S ratio were determined on samples 
taken from the same sheet; the modulus of elasticity varied be- 
tween 748,000 and 775,000 kg/cm 2 (average: 760,000 kg/cm 2 ); 
Poisson*s ratio, m = 3. 33. 

The wrinkling tests were made with the aforementioned clamp- 
ing device (Fig. 10). The effect of the edge loading on the dis- 
placement is graphed in Figure 32 (plate test No. 66). Tents 
42, 44, 45, 48, 52-59, and. 64-67 were made with the 0.03 cm 
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plates, and. the load at wrinkling and the wave length determined. 
The data are shown in Table VII and Figure 33. The graph shows 
the transition from elastic to inelastic deflection "between 
a/h = 53.3 and 33.3, and in conformity the k 1 values up to 
a/h =53.5 show a constant value of about 3,490,000. The fig- 
ures for k range between 51.1 and 47.6; the average is 49.2. 
Below a/h = 53.3 the .k l values continue to decrease to 
1,940,000. 

The experiments (Hos. 43, 46, 49-51, and 60-63) with the 
0 . 05-centimeter plates are included in Table VT.I, and Figure 34 
denotes the effect of a/h on the collapsing load. The transi- 
tion from elastic to inelastic deflection occurs between a/h = 
56 and 40.. The. k 1 values are a/;ain praotically constant, 
ranging between 3,740,000 and 3,500,000.. k averaged 51.0, a 
figure only slightly different frpm the. previously obtained 
value . .... 

Figure 35 .exhibits the \/a values for tests 42-67, plotted 
against plate width a. The a/h. ratio fluctuated between 
116.8 and .33.3, and the ratio of wave length to plate width be- 
tween (1.68) 1.89 and 2.02. The average for \/a for all tests 
with duralumin plates was \/a = 1.97. 
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TABLE VII. Tests 42-67 



Hat er i al , Dure 1'imi a 

Modulus of elasticity, E = 760,000 kg/cm a 

Poisson's ratio, m = 3.23 





Dimensions 




Wrink- 








ITT . _ 

Wave 




Test 


width 


thick- 


a/h 


ling 








length 


X 


No. 


a 


ness 
h 


load 
s 

kg/cm 




k 1 


K 


X 


a 




cm 












cm 




42 


2 . 00 


0 . 03 


66 . 6 


24.6 


3 , 


620,000 


50.9 


3 . 92 


1.96 


4-3 


2.00 


0 . 05 


40 


93 


2, 


676,000 


( 41 . 8 ) 


3 .96 


1 .98 


44 


1.50 


0.03 


53 3 


35.8 


3 . 


394,000 


47. 7 


3.17 


1.98 




2 . 80 


0 .03 


93 . 3 


12.1 


3^ 


513,000 


49. 4 


5.56 


1.985 


*±o 


2 . 80 


0 .05 


5G 


56 . 4 


3 . 


537,500 


49.8 


5.75 


2.05 


4R 


2 . 80 


0.03 


93 . 3 


12.2 


s, 


543,000 


49.9 


5.55 


1.98 


49 


2. 00 


0 .05 


40 


96 


3. 


072,000 


(43.2) 


3.95 


1.975 


50 


1 . 60 


0. 05 


33 


98 


2j 


007,000- 


(28.3) 


3.18 


1.985 


51 


1.60 


0. 05 


32 


103 


8, 


,110,000 


29.7 


3.10 


1.940 


52 


2.80 


0.03 


93 .3 


11.8 


3; 


,426,000 


48.2 


5.50 


1.960 


53 


2.80 


0.03 


93.3 


12.5 


3' 


,360,000 


51.1 


5.45 


1.945 


54 


2.00 


0.03 


bb . b 


OA A 


o . 


^ OJLO , uuu 


UU • o 


3 9,0 




55 


2.00 


0.03 


66.6 


24.3 


rr 

0 


,600,000 


50.6 


4.05 


2.02 


56 


1.60 


0.03 


53.3 


36. 7 


3' 


,480,000 


49.0 


3.10 


1.94 


57 


1.60 


0.03 


53,3 


35.7 


3. 


,385,000 


47.6 


3.18 


1.99 


58 


1.00 


0.03 


33.3 


52.5 


li 


,940,000 


(27.3) 


2.00 


2.00 


59 


1.00 


0.03 


33.3 


59 . 5 


2' 


,204,000 


(31.1) 


1.89 


1.89 


60 


2.80 


0.05 


56 


57.5 


1-7 
«J 


,746,000 


52.7 


5.39 


1.91 


61 


2.80 


0.05 


56 


55.5 


3' 


, 501 , 000 


49.6 


5.47 


1.95 


62 


3.50 


0.05 


70 


37 


3 


,626,000 


51.0 


7.15 


1.98 


63 


3. 50 


0.05 


70 


39 


3 


,822,000 


53.7 


7.05 


1.96 


64 


3.50 


0.03 


116 .8 


7.8 


3 


,547,000 


49.9 


7.20 


2.00 


65 


3.50 


0.03 


116.8 


7.2 


3 


,267,000 


46.0 


6.83 


1.95 


66 


1.60 


0.03 


53 . 3 


40 


3 


,792,000 


53.3 


2.79 


1.68 


67 


1.00 


0.03 


33.3 


55 


2 


,036,000 


(28.6) 


1.98 


1.98 
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TABLE VIII. Tests 68-80, 87, 88 

I^c^ "fc €/ 3? i G 1 1 j? C S 

Modulus of elasticity, E = 343,000 kg/cm 2 

Poisson's ratio, in = 3.06 

Plate thickness, h = 0-02 cm 



Test 
I>IO . 


L) "1 o "f" q 

r j.a o e 

width 
a 

Oiii 


a/ n 


W7I A V*l I/" — 

ling 
.Load. 

0 

Trier / pta 




y 
i- 


Wave 

"1 PTlCrfVl 
-L C 1 lc< Oil 

A. 

cm 


X. 
a 


68 


2.0 


100 


8.45 


4,225,000 


48.0 


4.35 


2.15 


69 


2.0 


100 


8.26 


4,130,000 


47.0 


3.80 


1.9 


70 


2.8 


140 


4.54 


4,449,000 


50.6 


5.46 


1.95 


71 


2.8 


140 


4.54 


4,449,000 


50.6 


5.43 


1.95 


72 


1.6 


80 


13.7 


4,334,000 


49.9 


3.17 


1.98 


73 


1.6 


80 


13.5 


4,320,000 


49.1 


3.38 


2.11 


74 


0.8 


40 


31.8 


2,544,000 


28.9 


1.54 


1.93 


75 


1.0 


50 








1.97 


1.97 


76 


1.2 


60 








2.25 


1.88 


77 


2.0 


100 


8.45 


4,225,000 


48.1 


4.08 


2.05 


78 


2.8 


140 


4.43 


4,341,000 


49.4 


5.38 


1.93 


79 


1.0 


50 


23.0 


2,875,000 


32.7 


1.90 


1.90 


80 


1.0 


50 


22.9 


2,863,000 


32.6 


1.82 


1.82 


87 


3.5 


175 


3.16 


4,933,550 


56.1 


7.52 


2.14 


88 


3.5 


175 


2.96 


4,625,000 


52 . 6 


6.3 


1.8 



Experiments with Brass Plates 



The samples were 0.02 cm thick rolled, brass. Tensile tests 
made on six specimens yielded E = 943,000 kg/cm 2 and transverse 
elongation m = 3.06. 

The results of the deflection tests are given in Table VIII 
and Figure 36. As long as ' k ' remains constant by a/h = 70, 
the deflection is elastic (tests 68-73, 77, 78, 87, 88). When 
a/h < 70, formula. (23) is inapplicable to this material; k» 
becomes smaller, so that hy a/h =40 it has dropped to 
2,544,000. For elastic deflection k = 48.0 - 56.1; the average 
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is k = 50.1. 

The ratio of wave length to plate width X/a, remains prac 
tically constant and raises between 1.8 and 2.15. The average 
of all tests is 1.96. In Figure 37, where \/a is plotted 
against the plate width, the individual X/a figures are in 
close agreement with the mean values shown as dotted line. 

Conclusion 

Owing to the high elastic def ormability of celluloid, it 
was not only possible to observe the beginning but also to ascer 
tain the type of deflection. The test data on celluloid were 
affirmed by the experiments with duralumin and brass. 

According to our experiments the load at wrinkling can be 

accurately, expressed by the formula s = k (k = 49 'to 50), 

a 

as set up from experiments with celluloid and substantiated by 
those on duralumin and brass. Its limit of applicability is 
the proportional elastic limit. 

In wrinkling, two cases must be distinguished: deflection 
prior to reaching yield limit (elastic deflection); and deflec- 
tion after the yield limit has been reached (inelastic deflec- 
tion). The transition from elastic to inelastic deflection in 
the individual materials is characterized by a certain ratio of 
plate width to plate thickness a/h, which for celluloid is 34, 
for duralumin, 50, and for brass, 70. The greater the elastic 
deflection in the material, the lower the a/h ratio of transi- 
tion. In contrast to Southwell and Skan, our figures for the 
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wrinkling loads were about 43$ lower. This discrepancy finally 
becomes greater than we expected from the previously "stated ef- 
fect of inexact edge constraint. Then the shape of the deflec- 
tion surely has some effect also. According to Southwell, the 
wave lengths should equal 1.6 times the plate width, but our 
tests yielded 1.97 a in every case. Timoschenko T s figures are 
slightly higher, although it should be remembered that in his 
tests all four edges v/ere supported. 

Regarding the- deflection itself, the following may be stat- 
ed: Once a certain load - the wrinkling load - has been ex- 
ceeded, the equilibrium between the internal and external 
stresses ceases to be certain. The beginning of instability is 
accompanied by a formation of regular, wavelike wrinkles. Near 
the free edge the uniformity is interrupted by the uneven 
stresses prevailing at that point. The deflection at the free 
edge is approximately 15$ greater than in the plate center. 
The wave length. of the deflection is in a constant ratio to the 
plate width, ; X/a ■= 1.97, and is unaffected by the plate thick- 
ness and the material. 

Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. ... 
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Fig. 25 Wrinkling by, s = 4.83 kg/cm, edge loading. 
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Fig. 28 Poriasgaont deformation in 0.05X1 . 79 ®m plate- 
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Figs. 29,30,31 
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Figs. 32, 33 
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Figs. 36,37 
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